Partial Wave Analyses of J/?/? ^ jK^K and 'jK^K^ 

J. Z. Bail, Y. Ban^o, J. G. Bian^, D. V. Bugg^^, X. Cai^, J. F. Chang^, H. F. Chen^'^ , H. S. Chen^ 
H. X. Chen^, Jie Chen^ J. C. Chen\ Y. B. Chen\ S. P. Chi\ Y. P. Chu\ X. Z. Cm\ H. L. Dai^ 
Y. S. Dail^ Y. M. Daf , L. Y. Dong\ S. X. Dul^ Z. Z. Du\ J. Fang\ S. S. Fang\ C. D. Fu\ 

H. Y. Fu\ L. P. Fu^ C. S. Gao\ M. L. Gao\ Y. N. Gaol^ M. Y. Gong\ W. X. Gong\ 
S. D. Gu\ Y. N. Guo\ Y. Q. Guo\ Z. J. Guo^^, S. W. Han\ F. A. HaW^ J. Re\ K. L. He\ 
g '. M. He^^ X. He^ Y. K. Heng\ H. M. Hu\ T. Hu\ G. S. Huang\ L. Huang^ X. P. Huang\ 

^ ; X. B. Ji\ Q. Y. Jia^o, C. H. Jiang\ X. S. Jiang\ D. P. Jin^, S. Jm\ Y. Jin\ Z. J. Ke\ 

Y. F. Lai\ F. Li^, G. Li^, H. H. Li^, J. Li^, J. C. Li^, K. Li", Q. J. Li^, R. B. Li^, R. Y. Li^ 
W. U\ W. G. Li\ X. Q. Li^, X. S. Li^^, Y. F. Liangl^ H. B. Liao^ C. X. Uu\ Fang Liu^^ 
F. Liu^ H. M. Liu\ J. B. Liu\ J. P. Liu^^ R. G. Lm\ Y. Liu\ Z. A. Liu\ Z. X. Liu\ G. R. Lu^ 



(N 
(N 



> ■ F. Lu\ H. J. Lu^^ J. G. Lu\ C. L. Luo®, X. L. Luo\ E. C. Ma^ F. C. Ma^ J. M. Ma\ 
00 i 

in ■ L. L. Ma^^^ X. Y. Ma\ Z. P. Mao\ X. C. Meng\ X. H. Mo\ J. Nie\ Z. D. Nie\ S. L. 01senl^ 
O 

^ . H. P. Pengi^ N. D. Qi\ C. D. Qian^^, J. F. Qiu\ G. Rong\ D. L. Shen^, H. Shen\ X. Y. Shen^ 

g : H. Y. Shengi, F. Shi\ L. W. Song\ H. S. Sun\ S. S. Sun^^ Y. Z. Sun\ Z. J. Sun\ S. Q. Tang^ 

■ X. Tangi, D. Tian^ Y. R. Tian^s, G. L. Tongi, G. S. Varneri^, J. Z. Wang^, L. Wang\ 
CD ■ 

ji^. L. S. Wang\ M. Wang^ Meng Wang^ P. Wang\ P. L. Wang\ W. F. Wang^ Y. F. Wang\ 

^ ; Zhe Wangi, Z. Wang^, Zheng Wang^, Z. Y. Wang^, C. L. Wei\ N. Wu^, X. M. Xia^, X. X. Xie^, 



X 



G. F. Xu\ Y. Xu\ S. T. Xue\ M. L. Yan^^ W. B. Yan\ F. Yang^ G. A. Yang\ H. X. Yangl^ 
J. Yangi^ S. D. Yang\ Y. X. Yang^, M. H. Ye^, Y. X. Ye^^ J. Ying^o, C. S. Yu\ G. W. Yu^ 
C. Z. Yuan\ J. M. Yuan^, Y. Yuan\ Q. Yue^ S. L. Zang\ Y. Zeng^ B. X. Zhang^, 
B. Y. Zhang\ C. C. Zhang\ D. H. Zhang\ H. Y. Zhang\ J. Zhang\ J. M. Zhang^, 
J. W. Zhangi, L. S. Zhang^, Q. J. Zhang^ S. Q. Zhang\ X. Y. Zhang^^, Yiyun Zhang^^ 
Y. J. Zhangio, Y. Y. Zhang\ Z. P. Zhang^^, D. X. Zhao^ Jiawei Zhao^^, J. B. Zhao^ 
J. W. Zhaoi, P. P. Zhao^, W. R. Zhao^, Y. B. Zhaoi, Z. G. Zhao^*, J. P. Zheng^, L. S. Zheng^, 
Z. P. Zhengi, X. C. Zhong^, B. Q. Zhou\ G. M. Zhou\ L. Zhou\ N. F. Zhou\ K. J. Zhu\ 
Q. M. Zhu\ Yingchun Zhu\ Y. C. Zhu\ Y. S. Zhu\ Z. A. Zhu\ B. A. Zhuang\ B. S. Zou^ 

(BES Collaboration) 
1 Institute of High Energy Physics, 
Beijing 100039, People's Republic of China 



1 



^ China Center of Advanced Science and Technology, 
Beijing 100080, People's Republic of China 
^ Guangxi Normal University, Guilin 54IOO4, People's Republic of China 
^ Henan Normal University, Xinxiang 453002, People's Republic of China 
^ Huazhong Normal University, Wuhan 430079, People's Republic of China 
^ Hunan University, Changsha 410082, People's Republic of China 
^ Liaoning University, Shenyang 110036, People's Republic of China 
^ Nanjing Normal University, Nanjing 210097, People's Republic of China 
^ Nankai University, Tianjin 300071, People's Republic of China 
Peking University, Beijing 100871, People's Republic of China 

Queen Mary, London EI4NS, UK 
Shandong University, Jinan 250100, People's Republic of China 
Shanghai Jiaotong University, Shanghai 200030, People's Republic of China 
Sichuan University, Chengdu 610064, People's Republic of China 
Tsinghua University, Beijing IOOO84, People's Republic of China 
University of Hawaii, Honolulu, Hawaii 96822 
University of Science and Technology of China, 
Hefei 230026, People's Republic of China 
Wuhan University, Wuhan 430072, People's Republic of China 
Zhejiang University, Hangzhou 310028, People's Republic of China 
* Visiting professor at the University of Michigan, Ann Arbor, MI 48109 USA 

(Dated: July 21, 2003) 



2 



Abstract 

Results are presented on J/^ radiative decays to K'^K~ and KgKg based on a sample of 
58M J /if; events taken with the BES II detector. A partial wave analysis is carried out using the 
relativistic covariant tensor amplitude method in the 1-2 GcV mass range. There is conspicuous 
production due to the /2(1525) and /o(1710). The latter peaks at a mass of 1740±4l25 MeV with 
a width of 166^gjl];Q MeV. Spin is strongly preferred over spin 2. For the /2(1525), the helicity 
amplitude ratios are determined to be = 1.00 ± 0.281^;^^ and = 0.44 ± O.OSlg;^^. 

PACS numbers: 14.40.Cs, 12.39.Mk, 13.25.Jx, 13.40.Hq 
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I. INTRODUCTION 



QCD predicts the existence of glueballs, the bound states of gluons, and the observation 
of glueballs is, to some extent, a direct test of QCD. Such gluonic states are expected to give 
rise to a rich isoscalar meson spectroscopy, and Lattice Gauge Theory calculations predict, 
in particular, that the lowest-lying state should occur in the mass range 1.4-1.8 GeV and 
have J^*" = 0"^"^ [l] . For a, J /ip radiative decay to two pseudoscalar mesons, only J^*^ values 
in the series O"*"^, 2"*"+, ... are possible, so such states provide a very clean laboratory to 
search for the lowest mass scalar glueball. 

There has been a long history of uncertainty about the properties of the /o(1710), one of 
the earliest glueball candidates. This history is reviewed in detail in the latest issue of the 
Particle Data Group (PDG) B] and will not be repeated here. The latest analysis of Mark 
III data by Dunwoodie B favors = 0^ over arr earlrer ass™,rt of 2^ while the latest 
central production data of WA76 and WA102 also favor O"*" j^|5|. In this paper, we present 
new results on J/ip —>■ '-/K^K^ and '-/KgKg based on a sample of 58M J/i/j events taken 
with the upgraded Beijing Spectrometer (BESII) located at the Beijing Electron Positron 
Colhder (BEPC). 



II. BES DETECTOR 

BES II is a large solid-angle magnetic spectrometer that is described in detail in Ref . . 
Charged particle momenta are determined with a resolution of (Jp/p = 1.78%-y/l +p2(GeV^) 
in a 40-layer cylindrical drift chamber. Particle identification is accomplished by specific ion- 
ization [dE/dx) measurements in the drift chamber and time-of-flight (TOF) measurements 
in a barrel-like array of 48 scintillation counters. The dE/dx resolution is (JdE/dx = 8.0%; the 
TOF resolution is axoF = 180 ps for Bhabha events. Outside of the time-of-flight counters 
is a 12-radiation-length barrel shower counter (BSC) comprised of gas proportional tubes 
interleaved with lead sheets. The BSC measures the energies and directions of photons with 
resolutions of <7e/E ~ 21%/^i?(GeV), cr^ = 7.9 mrad, and o"^ = 2.3 cm. The iron flux 
return of the magnet is instrumented with three double layers of counters that are used to 
identify muons. The average luminosity of the BEPC accelerator is 4.0 x 10'^° cm^^s^^ at 
the center-of-mass energy of 3.1 GeV. 
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In this analysis, a GEANT3 based Monte Carlo simulation package (SIMBES) with de- 
tailed consideration of real detector performance (such as dead electronic channels) is used. 
The consistency between data and Monte Carlo has been carefully checked in many high 
purity physics channels, and the agreement is quite reasonable. 



III. EVENT SELECTION 



The first level of event selection requires two charged tracks with total charge zero for 
'~fK~^K~ candidate events, and requires two positively-charged and two negatively- charged 
tracks for 'jK^Kg events. These tracks are required to lie well within the acceptance of the 
detector and to have a good helix fit. More than one photon per event is allowed because 
of the possibility of fake photons coming from the interactions of charged tracks with the 
shower counter or from electronic noise in the shower counter. 

For J/ip ^ '~fK'^K~, the vertex is required to lie within 2 cm of the beam axis {x — y 
plane) and within 20 cm of the center of the interaction region (along z). Each of the 
charged particles is required to not register hits in the muon counters in order to remove 
7/i"^/i~ events. The following selection criteria are used to remove the large backgrounds 
from Bhabha events: (i) The opening angle of the two tracks satisfies 9 op < 175°. (ii) 
The energy deposit of each track in BSC satisfies Esc < 1-0 GeV. In order to reduce the 
background from final states with pions and electrons, each event is required to have at least 
one kaon identified by the TOF. Requirements on two variables, U and P^^, are imposed 
0]. A "missing- neutral-energy" variable U = {Emiss — I Pmiss I) is required to satisfy 
—0.10 < U < 0.20 GeV; here E^iss and Pmiss are the missing energy and momentum of 
all charged particles respectively. Also a "missing-p^" variable P'l^ = 4| Pmiss psin^ 6*^/2 is 
required to be < 0.002 GeV^, where 6^ is the angle between the missing momentum and 
the photon direction. The U cut removes most background from events having multipion 
or other neutral particles, such as pn, ■yn~^n~ events; P^^ is used to eliminate background 
photons. The selection criteria for a good photon used here are based on those applied in 
previous BESI analyses 8]. In brief, the good photon is required to be isolated from the 
two charged tracks and to come from the interaction point. 

In order to reduce the J/ip —>■ tt^K^K^ and J/ip —>■ tt'^tt+tt^ contamination, all events 
surviving the above criteria which have two or more photons are kinematically fitted to 
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M(KK) (GeV/c^) 

FIG. 1: Invariant mass spectra of a) K^K^ , b) K'^Kg for J/ip —>■ ^KK events, where the shaded 
histograms correspond to the estimated background contributions. 

these hypotheses. Those events with a fit < 50, and with photon pair invariant mass 
within 50 MeV/c^ of the 7r° mass, are rejected. Finally, the two charged tracks and photon 
in the event are 4-C kinematically fitted to obtain better mass resolution and to suppress 
backgrounds further by the requirements x^k+k- < ^^"^ X^k+k- < 'X^-w+-k-- 

For J/ip ^ 'jKgKg, the Kg mesons in the event are identified through the decay Kg — >■ 
TT+TT^. The four charged tracks can be grouped into two pairs, each having two oppositely 
charged tracks with an acceptable distance of closest approach. Signal events are required 
to satisfy 6]^^ < (20MeV/c2)2, where 6]^^ = (M^+^-(l) - MksY + (M^+^-(2) - Mks? and 
M^+^- is calculated at the Kg decay vertex. The main backgrounds from '-fKgK^n'^ and 
"yKgKgii^ events are suppressed by requiring [/ < 0.10 GeV, P^^ < 0.005 GeV^ and the 4-C 
kinematic fit x^in < 10- 

Fig. ^ shows the K^K^ and KgKg mass spectra for the selected events, together with 
the corresponding background distributions. These two mass spectra agree closely below 
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2.0 GeV. The resonant structures in the mass regions of the /2(1525) and the /o(1710) 
are very clearly visible in both decay modes. Averaged over the whole mass range, the 
detection efficiency for 'yK~^K~ is 14.7% and for 'jKgKg is 14.5%. For the 'yK~^K~ channel, 
the experimental background arises mainly from the non-resonant K^K^n^ and two-body 
K*^K^ events which are peaked at high K^K^ masses. In the entire mass range, 14597 
^K^K^ events are reconstructed, and the detailed Monte Carlo simulation of the BES 
detector estimates a background of 3094 events. The estimation of the background events in 
the ^K^Kl sample is obtained from the side band (28.7MeV/c2)2 < < {?,bMeV/c^f- 
this equal-area-selection provides a properly normalized background estimation. In Fig.^), 
there are 3169 selected -^K^K^ events and 413 background events. 



IV. ANALYSIS RESULTS 



We have carried out partial wave analyses using relativistic covariant tensor amplitudes 
constructed from Lorentz-invariant combinations of the 4-vectors and the photon polariza- 
tion for J/ip initial states with helicity ±1 [^. Cross sections are summed over photon 
polarizations. The relative magnitudes and phases of the amplitudes are determined by a 
maximum likelihood fit. The background events obtained from Monte Carlo simulation or 
5\-g side band are included into the data samples, but with the opposite sign of log likelihood 
compared to data. These events cancel background within the data samples. These analyses 
are confined to masses less than 2 GeV in order to ensure that a description containing only 
0^^ and 2++ amplitudes be appropriate. The mass distributions of K^K~ and K^K^ after 
acceptance and isospin corrections for missing '-yK^K^ and ""jK^K^ with -k^ti^ decays 

are shown in Fig. |21 The event topologies of the K^K~ and K^K^ modes are different, so 
that acceptance and background effects are rather different also; nevertheless, there is good 
quantitative agreement between the two distributions. 



A. Bin-by-bin analysis 

In the bin-by-bin analysis, the data in mass intervals 40 MeV wide are ffited with four 
helicity amplitudes, one for = 0"*" and three for 2^ amplitudes Q]. The mass interval 
width is chosen as a compromise between the desire for high statistics in each mass interval. 
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FIG. 2: The mass distributions of KK from J/ip radiative decays, after acceptance and isospin 
corrections for missing ^K^K^ and ^K^K^ with vr^vr'^ decays. 

and the need for detailed information on the mass dependence of each measured amphtude. 
In each mass interval, the jKR data sample is analyzed in terms of the joint production and 
decay angular distribution of the pseudoscalar meson system. The S- and D-wave intensity 
distributions, |ao,oP; k2,oP) ^2,1^ and |a2,2p for jKK data resulting from this bin-by-bin 
fit are shown function of mass in Fig. El 

The KK S-wave intensity dominates the 1.7 GeV region. The solid curves in Fig. El 
correspond to fits of coherent superpositions of individual Breit-Wigner resonances to the 
data points of each intensity distribution. The following channels are considered: 



7/2(1525) 
7/o(1710) 
7/2(1270) 
7/o(1500) 
7 + broad 0" 



and 2^^ components 



The first two are dominant. There is evidence for existence of the /2(1270), and the /o(1500) 



10000 




M(KK) (GeV) 

FIG. 3: The mass dependence of the amphtude intensities for ^KK data. The sohd curves 
correspond to the coherent superposition of the Breit-Wigner resonances fitted to the acceptance- 
and isospin-corrected data points obtained from the bin- by-bin fit. The dashed Hne histograms are 
the results of the global fit described in the text. 



is included here for consistency with the global fit below. 

For the spin amplitude, two interfering resonances (/o(1500, /o(1710)) and an interfering 
constant amplitude term, which is used to describe the broad S- wave contribution, are 
included. The mass and width of the /o(1500) are fixed to the PDG values; those of the 
/o(1710) are to be determined. The /o(1710) is well described by a Breit-Wigner of mass and 
width M = 1722 ± 17 MeV, T = 1671^^ MeV, and the branching fraction for J/i) radiative 
decay to the combined KK modes is i3(J/^ 7/o(1710) 'yKK) = (ILl^La) x 10"^. 
The errors here are statistical errors. 

For the spin 2 amplitudes, the /2(1525) and /2(1270) are included. There is also some 2++ 
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structure above 2.0 GeV in KK mass, which could contribute to the present fitted range, 
and thus the tail of a high mass 2++ state is included in our fit. We choose a resonance mass 
of 2250 MeV and width of 350 MeV to represent the structure in the higher mass region. The 
mass and width of the /2(1270) are fixed at the values quoted in the PDG. For the tensor 
resonance, /2(1525), its mass and width are fixed to the values M = 1519 MeV, F = 75 
MeV determined by the global fit which is described below, and the total branching fraction 
and ratios of amplitude intensities are determined to be B{J/ip — > 7/2(1525) jKR) = 
(4.02 ±0.51) X 10"^; = |a2,i| Vl«2,or = 1-32 ±0.29, = |a2,2| Vl«2,or = 0.38 ±0.20. The 
intensity of the /2(1270) is poorly measured because of the relatively low statistics and the 
weak coupling of this state to KK. The amount of spin 2 component in the 1.7 GeV mass 
region is small, ~ (16 ± 9)%. The errors shown above are statistical and are obtained from 
the Breit-Wigner fit. 

B. Global fit analysis 

We now turn to the global fit to the J/ip ^ 'yK~^K~ and J/ip ^ 'yK^K^ data. Each 
sample is analyzed independently, and the fit results shown below are for their averaged 
values. This fit has the merit of constraining phase variations as a function of mass to simple 
Breit-Wigner forms. It also performs the optimum averaging of helicity amplitudes and 
their phases over resonances. Partial waves are fitted to the data for the same components 
described in the bin- by-bin fit. The broad 0"*""*" component improves the fit significantly; 
removing it causes the log likelihood value to become worse by 221. For the /2(1270) and 
/o(1500), we use PDG values of masses and widths, but allow the amplitudes to vary in the 
fit. For the /2(1525), relative phases are consistent with zero within experimental errors. It 
is expected theoretically that relative phases should be very small, on order of a ~ 1/137 for 
the electromagnetic transitions J/ip ^ 'y + 2"*". In view of the agreement with expectation, 
these relative phases are set to zero in the final fit, so as to constrain intensities further. 

A free fit to /2(1525) gives a fitted mass of 1519 ± 2 MeV and a width of 75 ± 4 MeV. 
The fitted mass and width of the /o(1710) are M = 1740 ± 4 MeV and F = 166^8 MeV, 
respectively. The fitted intensities are illustrated in Fig. HI For the /2(1525), we find the 
ratios of helicity amplitudes x"^ = 1.00 ±0.28 and y"^ = 0.44 ±0.08. In this fit, we allow some 
O"*" contribution under the /2(1525) peak, while previous analyses by DM2 and Mark III 
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FIG. 4: The KK invariant mass distributions from J/i/j ^K^K~ and J/ifj ^ ^K'gK'g. The 
points are the data and the full histograms in the top panels show the maximum likelihood fit. 
Histograms on subsequent panels show the complete 0"*" and 2+ contributions including all inter- 
ferences. 

[10, 11] ignored the small 0"^ contributions. The branching fractions of the /2(1525) and the 
/o(1710) determined by the global fit arei3(J/V^ ^ 7/2(1525) ^ -fKK) = (3.42±0.15) x 10"^ 
and B{J/^ 7/o(1710) -fKK) = (9.62 ± 0.29) x 10"^ respectively. The errors shown 
here are also statistical. An alternative fit to /j(1710) with = 2+ is worse by 258 in log 
likelihood relative to O"*" for '-/K^^K" data and by 67 for '-/K^Kg. Remembering that three 
helicity amplitudes are fitted for spin 2 but only one for spin 0, the fit with = 0"^ is 
preferred by > lOa after considering the two data samples together. 

The separation between spin and 2 is illustrated in Fig. taking the J/ip ^ 'yK~^K~ 
data as the example. Let us denote the polar angle of the kaon in the KK rest frame by 
6k, and the polar angle of the photon in the J/ip rest frame by 6-y. The data are fitted 
simultaneously including important correlations between 6k and 6^. The left panels show 
resulting fits to cos 6k for J = and 2. There is no significant difference between the two 
fits. The distributions should be flat for O"*", but the interference with the tail of /2(1525) 
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FIG. 5: Projections in cos 9k and cos^-y for 0^^ and 2^^ assumptions. The points are the data 
{J/tl' — > jK^K^ sample), and the histograms are the global fit results. 

has a large effect. The right panels show the fits to cos 6'^; the optimum fit is visibly better 
for J = than for J = 2. [If one fits only the cos^^ distribution, it is possible to fit equally 
well with J = or 2, but then the fit to cosOk gets much worse.] 

If the /o(1500) is removed from the fit, the log likelihood is worse by 1.65 (3.58) for K^K~ 
{K'gK'g)^ corresponding to about 1.3cr (2.2cr). If the /2(1270) is removed, the likelihood is 
worse by 57.5 (13.6) for K^K^ {KgKg), corresponding to > 5cr (3.8(t). 

V. SYSTEMATIC ERROR 



The systematic error for the global fit is estimated by adding or removing small compo- 
nents used in the fit, replacing the /o(1500) with the /o(1429), V = 169 MeV, described in 
Ref. P|, varying the mass and width of the large /2(1525) within the PDG errors, varying 
the mass and width of /o(1710) based on the difference between the K^K~ and K^K^ decay 
modes, and varying the background component within reasonable limits in both the global 
fit and bin-by-bin fit. It also includes the uncertainty in the number of J /ip events analyzed 
and the difference from two different choices of MDC wire resolution simulation. 

The uncertainty about the shape of broad O"*""^ background is included in the systematic 
error also. An incoherent fit with this broad component and a fit with alternative forms for 

n 

the s-dependence using the parametrization of Zou and Bugg II4I for the /o(400 — 1200) 
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have been performed to estimate the systematic error from this source. This uncertainty 
affects the results significantly, especially the branching fractions, because of the interference 
between the broad structure and the other components. Therefore, the error from this 
model-dependence for the branching fraction measurements is separated from the statistical 
and other systematic errors in our final results. The systematic errors for the global fit are 
summarized in Table HI For the mass and width, only the contributions from the model- 
dependence, which are large compared to the other errors, are shown in the table. 





M/^(1525) 


r/^(1525) 






^/'(1525) ^MniO) r J (1710) 


^/o(1710) 


remove /o(1500) 






+32 
—0 


+20 
—0 


±0 


+10 
—0 


use /o(1429) 






+0 
-15 


+0 
-9 


+0 
-5 


+3 

-0 


remove /2(1270) 






+42 
-0 


+0 
-55 


+6 
-0 


+0 

-1 


use the a 


+0.66 


+20 


+17 

-9 


+0 
-14 


tf -1.44 +9 


+29 
-0 


incoherent 0^^ 


+0.99 




+6 
-0 


+0 
-64 


tf - +3 


+28 
-0 


M, r of /^(1525) 






+49 
-15 


+0 
-34 


+11 

-8 


+4 
-5 


M, rof /o(1710) 






+51 
-17 


+11 

-36 


±3 


+1 
-0 


M, r of high 2++ 






+46 
-14 


+0 
-59 


+1 

-4 


+6 
-0 


background 






+46 
-17 


+0 
-55 


+0 
-3 


+9 
-10 












±4.7 


±4.7 


wire resolution 










±15 


±15 



TABLE I: Estimation of systematic error (%) in the global fit. -6^^(1525) and Sj(,(i7io) are the 
branching fractions for /2(1525) and /o(1710) respectively. 



VI. RESULTS AND DISCUSSION 

The results of the bin-by-bin and global fits are summarized in Tables |n] and respec- 
tively. For the bin-by-bin fit, the errors are statistical ones only, and for the global fit, the 
first error listed is the statistical error, the second error is the systematic error, and the third 
one for the branching fractions is for the model-dependence of the broad components. 

The two fit methods, bin-by-bin and global, are based on different analysis concepts. 
In the bin-by-bin fit, the S- and D-wave intensities are fairly well determined and nearly 
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/^(1525) 


/o(1710) 


M (MeV) 


1519 (fixed) 


1722 ± 17 


r (MeV) 


75 (fixed) 




B( J/V' ^ iX, 


4.02 ±0.51 


'-^■^-1.2 


X = |a2,i! /|a2,ol 


1.32 ±0.29 




= |02,2|^/|02,OP 


0.38 ± 0.20 





TABLE II: Measurements of the /2(1525) and /o(1710) for the bin-by-bin fit. Errors shown are 
statistical only. 





/^(1525) /o(1710) 


M (MeV) 


1519 ± 2tl^ 1740 ± Atl^ 


r (MeV) 


75 ± 4l^^ 166i|ll^ 


X Jv J?)(xlO-^) 


■1 12 + n 1 r,+0.69+1.55 q , ^ .^n+2.11+2.Sl 
,). iZ i U- J--J-().(,5-0.()() ^ "■-•-'-1. 80-0. 00 


amp. ratios 


1.00 ± 0.28iJ:i]^ 




0.44 ± o.os^H:^^ — 



TABLE III: Measurements of the /2(1525) and /o(1710) for the global fit. The first error is 
statistical, the second is systematic, and the third is that corresponding to model-dependence of 
the broad components. 

model independent. The only model dependence in the bin-by- bin fit is the assumption 
that only S- and D-waves need be considered; this is reasonable, since one would not expect 
significant 4+"*" amplitudes below 2 GeV. However, due to limited statistics for each bin and 
the limited solid angle coverage of the detector, the relative phases of partial waves cannot 
be well determined. This causes larger uncertainties when extracting the mass and width of 
resonances by fitting only the partial wave intensities without the constraints of the relative 
phases between them. In the global fit, the phase variations as a function of mass are 
constrained to simple Breit-Wigner (BW) forms . The stability of the minimum optimizing 
procedure and statistical errors arc better than those of the bin-by-bin fit. However, if some 
non-BW resonance is assumed to be a BW-form amplitude, this will give a model-dependent 



14 



biased result. The model independent bin-by-bin result for the partial wave intensities can 
provide guidance for choosing components for the global fit. The final full amplitudes from 
the global fit definitely give a better fit to the whole set of data than the amplitudes obtained 
from fitting the partial wave intensities without constraints of relative phases between them. 

Fortunately from Tables [H] and 11111 and the comparison shown in Fig. |21 we see that 
the results obtained from the bin-by-bin fit and the global fit for the /2(1525) and /o(1710) 
agree with each other well within the errors. The ratios of the helicity amplitudes of the 

a 1525) from the present analysis are in reasonable agreement with Krammer's predictions 
. These ratios provide useful information for testing models of the resonance production 
and decay mechanisms. Most importantly, the analysis demonstrates that the mass region 
around 1.7 GeV is predominantly 0^^ from the /o(1710) Q]; this conclusion is consistent 
with that of references [3-5]. 



VII. SUMMARY 



In summary, the partial wave analyses of J/ip 'yK~^K~ and J/ip ^ 'jKgKg using 
58M J/ip events of BES 11 show strong production of the /2(1525) and the S-wave resonance 
/o(1710). This confirms earlier conclusions that the spin-parity of the /o(1710) is = 
0+. The /o(1710) peaks at a mass of 1740 ± 4^^^ MeV with a width of I66tlt\l MeV. 
For the /2(1525), the hehcity amplitude ratios are determined to be 1.00 ± 0.28to;36 and 
0.44 ± 0.08to;5g, respectively. They are consistent with theoretical predictions. 
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